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ABSTRACT 

We explore the elemental abundance features of metal-rich disk stars, highlighting the comparisons 
made with those of the recently revealed Galactic bulge stars. A similarity between two of the 
comparisons leads to a new theoretical picture of the bulge-disk connection in the Galaxy, where a 
supermassive black hole resides at the center. We postulate that a metal-rich outflow, triggered by 
feedback from a black hole, was generated and quenched the star formation, which had lasted several 
billion years in the bulge. The expelled gas cooled down in the Galactic halo without escaping from 
the gravitational potential of the Galaxy. The gas gradually started to accrete to the disk around 
five billion years ago, corresponding to the time of sun's birth, and replaced a low-metallicity halo gas 
that had been accreting over nearly ten billion years. The metal-rich infalling gas, whose elemental 
abundance reflects that of metal-rich bulge stars, mixed with the interstellar gas already present in the 
disk. Stars formed from the mixture compose the metal-rich stellar disk. This scheme is incorporated 
into models for chemical evolution of the disk. The resultant elemental features are compatible with 
the observed abundance trends of metal-rich disk stars, including the upturning feature exhibited 
in some [X/Fe] ratios, whose interpretation was theoretically puzzling. Furthermore, the predicted 
abundance distribution function of disk stars covers all observational facts to be considered: (i) the 
deficiency of metal-poor stars, (ii) the location of peak, and (iii) the extended metal-rich tail up to 
[Fe/H]~ +0A. 

Subject headings: Galaxy: bulge — Galaxy: center — Galaxy: disk — Galaxy: evolution — ISM: jets 
and outflows — stars: abundances 
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1. INTRODUCTION 

Feedback from supermassive black holes at the in- 
nermost regions of galaxies has a c rucial influence o n 
the evolution of their host galaxies (I Silk fc ReesI I1998D . 
While the growth of black holes through gas accretion 
is regulated by its associated feedback, the feedback 
also influences the host on a large scale through en- 
ergy input into the gas surrounding the black holes, 
which hal ts star formation an d induces a powerful out- 
flow (e.g.,l Springcl et al.' 2005) . Recent theoretical works 
(e.g. jPi M atteo ct al. 2005; Schawinski ct al. 2006) sug- 
gest that the processes related to feedback lead to the 
observed relationship between the black hole mass and 
the v elocity dispersion of th e galaxy spheroidal compo- 
nent (Magorrian et all Il998t iFerrarese fc MerrittI l2000l : 
iGebhardt ct al. 200^ 

It is evident that the Galaxy possesses a supermassive 
black h ole, with a rnass of ^ 3 x lO^Mp, at t he Galaxy 
center ()Genzel et al.lll997l: [Schddel et aLlbool . Because 
of the short distance, the Galaxy is the best-suited lab- 
oratory to probe the past of other galaxies hosting su- 
permassive black holes and that have been affected by 
their feedback. Now is the time to pursue this subject 
using a powerful tool to trace the mechanisms of galaxy 
formation, e.g., elemental abundances of stars, especially 
since we have obtained detailed elemental abundances of 
individual stars in the Galaxy, including member stars 
belonging to the bulge. 

The Galactic bulge is old (e.g.. lOrtolani et al]|1995D . 
No clear evidence of young stellar population can be de- 
tected in t he color-magnitude diagra m of the Galactic 
bulge stars (jFeltzing fc Gilmordl2000t) . although elemen- 



tal abundances of bulge stars (jFulbright et al.ir2007f ) with 
a signature of Type la supernovae (SNe la) imply that 
the period of star formation should be in the order of a 
billion years. In any event, the process of expelling the 
gas from the bulge to cease star formation at a young age 
of the Galaxy is required at the final stage of bulge for- 
mation, which can be identified with outflow energized 
by any past activity in the existing supermassive black 
hole. 

Outflow from the Galactic bulge should be well en- 
riched because its elemental feature is the end result 
of chemical evolution in the Galactic bulge. Thus, the 
metallicity of outflow will be almost equivalent to the 
highest among the bulge st ars, whose metaUicitie s range 
from -1.3 <[Fe/H]< +0.5 (iFulbright et al.l[200l . Such 
a metal-rich outflow crucially influences the chemical evo- 
lution of the Galaxy because of its high metallicity, as 
long as the gas released from the Galactic bulge does not 
have sufficient energy to escape from the gravitational 
potential of the Galaxy. This may be the likely fate of 
any outflow originating from a supermassive black hole 
with a relatively small mass like ours. 

Recently, detailed elemental ab undances of the Galac - 
tic bulge stars have been revealed (|Fulbright et al.ir2007f l. 
One of the most remarkable abundance features is the 
enhancement of [a-element/Fe] in comparison with disk 
stars in the overall metallicity range, including metal- 
rich stars above solar metallicity. These characteristics 
can be used for chemical tagging of an outflow. Consid- 
ering the cooling process that acted on the outflow while 
moving in the Galactic halo, a promising site where we 
could search for this tag would possibly be the Galactic 
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disk, although the process to gain angular momentum 
should be validated by future numerical simulations. 

Disk stars^ present a puzzling problem associated with 
their elemental feature for [Fe/H] 0. It is well un- 
derstood that the observed decreasing [a-element/Fe] 
trends with an increasing [Fe/H] for the range of 

— 1 [Fe/H] ;S are caused by a delayed Fe supply 
from SNe la compared with Type II SNe (SNe II). 
On the contrary, a theoretical interpretation of their 
abundance trends for [Fe/H] has not yet been ob- 
tained. Strangely, the [a-element/Fe] ratios, except for 
[0/Fe], start to upturn from [Fe/H]~ and continue to 
rise until reaching t he highest metallici ty [Fe/H]~ -1-0.4 
(|Reddv et all 120031 iBensbv et all l2005f ) . Such upturns 
can also be seen in other elements, such as odd ele- 
ments (Na, Al), and some iron-group elements (Ni, Zn). 
It is incorrect to attribute these upturns to metallicity- 
dependent SN II yields, especially because this feature 
is observed with both even- and odd-numberd elements. 
Alternatively, some new insights from the viewpoint of 
galaxy formation should be introduced to understand the 
late stage of chemical evolution in the Galactic disk im- 
printed in elemental abundances of metal-rich disk stars. 

Models of chemical evolution require that the Galac- 
tic disk is formed by a continuous infall of gas fro m the 
halo, to solve the so-called G-dwarf problem (e.g.. iPageH 
|1997[ ). Conventionally, it is assumed that an infalling 
gas has very low-metallicity over the entire period of disk 
formation. In practice, models adopting this assumption 
apparently reproduce the observed abundance distribu- 
tion function (ADF) of disk stars, as attempted by man y 
authors fe.g.. lYoshii et al.lll996t IChiappini et al.lll997D . 
However, as will be discussed in the next section, there 
exist a fatal problem with the present models when try- 
ing to reproduce the overall ADF by paying attention to 
not only its metal-poor tail but also its metal-rich side. 

Here, we first propose that there exists a switchover 
from accretion of low-metallicity halo gas to that of 
metal-rich gas onto the disk at the time when the metal- 
licity in the disk nearly reaches solar metallicity. This is 
followed by the scenario that a metal-rich gas blown away 
from the Galactic bulge cooled down radiatively in the 
halo and finally settled down on the disk that had been 
formed from a low-metallicity infall up until that time. 
Based on this theoretical scheme, the model of chemical 
evolution is designed and found to successfully reproduce 
the observed abundance trend of disk stars in the range 

— 1 <[Fe/H]< -1-0.4. Our scenario is supported by the 
presence of stars with super solar metallicities such as 
[Fe/H]= -1-0.2 — 0.4. We will begin with a discussion on 
this matter. 

2. IMPLICATIONS FROM ABUNDANCE DISTRIBUTION 
FUNCTION OF DISK STARS 

The distinct features of the ADF of solar neigh- 
borhood disk stars, confirmed by many studies thus 
far (e. g.. IWvse fc GilmOT3ll995t iRocha-Pinto fc Maciell 
Il996t iNordstrom et al.l l2004D are (i) the deficiency of 
metal-poor stars: there are few stars below [Fe/H] ;$ — 
1, i.e., the G-dwarf problem, and (ii) the peak lo- 
cated at [Fe/H]pcak"^ —0.2 — —0.1. In addition 

^ In this letter, disk stars include only stars belonging to the thin 
disk component in the Galaxy. 
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Fig. 1. — Abundance distribution function of disk stars against 
the iron abundance. Solid curves are the predictions of the model 
with [Fe/H]present = "1-0.4 and ti^=1.5 Gyr for tin=3 Gyr (left) 
and tin=5 Gyr (right), respectively. Results of the model with 
tla=1.5 Gyr and tin=5 Gyr for [Fe/HJprcsont = -1-0.03 and +0.15 are 
shown by dashed curve and dotted curve, respectively. These cal- 
culated distributions are convolved using the Gaussian with a dis- 
persion of 0. 15 dex in [Fe/H]. Filled and o pen circles represent dat a 
taken from Edvardsso n et al.l II1993I ') and IWvse fc Gilmor j II1995I ). 
respectively. The model distributions and the observed one by 
iWvse fc Gilmora il995i ) are normal ized to coincide with the total 
number of the sample stars used by lEdvardsson et ahl II1993I) . 



to these observational facts, spectroscopic observa- 
tions of elemen tal abundance s for metal-r ich disk stars 
(Feltz ing fc Gus tafsson 199^; iBensbv et a l. 2005) have 
revealed that chemical enrichment in the solar neigh- 
borhood has continuously proceeded until [Fe/H]~ -1-0.4. 
This information, relating to how far the stellar metallic- 
ity extends to a metal-rich direction, should be the third 
ingredient to be highlighted in the reproduction of the 
observed ADF. 

Next, we will try to reproduce the ADF of disk stars, 
paying special attention to the metal-rich side, which 
has not been fully considered by the previous studies. 
For this purpose, the procedure t o solve the equatio ns 
for chemical evolution adopted bv lYoshii et al.l (|1996f) is 
used, because in this approach, the present metallicity 
[Fc/HJtq of gas, which corresponds to the highest metal- 
licity of disk stars, is not an outcome but an input pa- 
rameter. It means that the value of [Fe/H]T'G can be 
chosen arbitrarily. This, together with the present gas 
fraction /^(Tg), is combined with four theoretical inputs 
composed of (i) the power index n of the star formation 
rate (SFR) in Schmidt law, (ii) the timescale tin of in- 
fall, (iii) the lifetime tu of SN la progenitors, and (iv) 
the metallicity [Fe/H]infaU of the infalling gas, to deter- 
mine the chemical evolution of the disk. For (i), it is 
assumed that the SFR is proportional to the gas den- 
sity, i.e., n = 1. For (iv), an infall from the halo under 
consideration implies a very low metallicity, likely to be, 
[Fe/H]i„faii^ —1 reflecting the metallicity of halo stars. 
Since this value is irrelev ant to the predictio n of ADF, as 
long as [Fe/H] infall < —1 (jYoshii et al.lll996[ ). we assume 
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[Fe/H]i,ifaii — 0. In the end, each combination of and 
tia yields various forms of the ADF. 

The sohd curves in Fig. 1 show the resultant 
ADFs with the choices of tia=1.5 Gyr and tin=3, 5 
Gyr, together with [Fe /H]rr.=+0.4 and /g(rG)=0.15 
() Young fc ScovilldfTOQlh . Here [Fe/H]TG is assumed to 
equal [Fe/H] of th e most metal-rich local F/G stars 
(|Bensbv et al.ll2005l ) . The age Tq of a galaxy is set to be 
14 Gyr. Note that we convolve the calculated ADF us- 
ing the Gaussian with a dispersion of 0.15 dex in [Fe/H], 
which mimics the error of the data. Obviously, the pre- 
dicted ADFs are biased to metal-rich as compared with 
the observations. Any combinations of the four input 
parameters, together with an allowable uncertainty in 
the observed /^(Tg), are unable to avoid positioning the 
peak beyond solar metallicity as long as the observed 
value for [Fe/H]TG is adopted. We performed another 
calculation with the assumption of [Fe/H]G=-I-0.15. This 
value corresponds to the mean metallicity at present 
in the observed age-metallicty relation of nearby stars 
(jNordstrom et al. 2004 : Reid et al. 2007). The resultant 
ADF is shown to be still shifted to metal-richer than the 
observation (dotted curve) . On the other hand, if we as- 
sume a smaller [Fb/HJtg, such as -1-0.03 (dashed curve), 
the calculated ADF matches the obs erved one , as already 
shown in the previous studies (e.g., [Yoshii_et al. 1996). 

In summary, the simultaneous reproduction of both 
[Fe/H]TG> -1-0.15 and [Fe/H]peak < is hard to realize 
through the conventional scheme adopted for chemical 
evolution of the Galactic disk. Thus, some modification 
of the constructed model is required. 

3. ABUNDANCE TRENDS OF METAL-RICH DISK STARS 

In this section, we will examine the abundance trends 
of metal- rich disk stars with [Fe/H] to persuade read- 
ers that the key factor required to solve the problem 
in the existing models is tightly linked to the activity 
and formation of the Galactic bulge. Figure 2 shows the 
comparison between abundance features in the [X/Fe] 
vs. [Fe/H] diagram for six elements of disk stars and 
metal-rich bulge stars. As already mentioned in Section 
1, the characteristic feature of a-elements, such as Mg, 
Si, Ca, and Ti, in comparison with Fe for metal-rich disk 
stars, is the upturn with an increase in [Fe/H]. Two up- 
per panels give a sample of X=Mg and Si, demonstrating 
that this feature is surely observed in ;$ [Fc/H] ;S +0.4. 
Although the origin of these upturn has been a totally 
open question, superposition of the observed data for 
metal- rich bulge stars on those of the disk stars leads to a 
possible indication for solving the puzzling upturn. Our 
proposed view of abundance trends between two differ- 
ent populations is that the [Mg (Si) /Fe] ratio of disk stars 
gradually approaches that of the metal-rich bulge stars 
in relation to the increase in [Fe/H]. This interpretation 
also holds true in the cases of [0/Fe] and [Mn/Fe], where 
the abundance trend in a metal-rich regime inherits that 
for [Fc/H] £ (middle panels of Fig. 2). In addition, 
odd elements such as Na and Al, which show an upturn 
in their ratios to Fe like the case of a-elements, also have 
enhanced ratios of bulge stars ahead in the direction of 
their upturn trends (lower panels of Fig. 2). 

However, this poses the question of which theoretical 
scheme can be applied to the potential bulge-disk connec- 
tion to understand their stellar elemental relationship. A 
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Fig. 2. — Correlations of [X/Fe] with [Fe/H] for disk stars (open 
circles and crosses) together with metal-rich ([Fe/H]> 0) bulge 
stars (filled circles). Except for Mn, filled and open circles rep- 
resent data from Fulbright et al. (2007) and Bensby ct al. {200^Y 
respectively. The Mn abundances are taken f r om|McWilliam et all 
(2003 ) (filled circl es), Feltzing fc Gustaf'ssonI l ll998l) (open circles'). 
and lReddv et af] (12003. ) (crosses). 



possible answer is that metal-rich disk stars are formed 
from a mixture of an infalling metal-rich gas originally 
ejected from the bulge and the remaining gas of a star 
formation in the disk. In this case, their abundance pat- 
tern is determined by the combination of heavy elements 
ejected from the bulge at the end of star formation, along 
with those elements that are already present in the in- 
terstellar gas of the disk. In accordance with an accumu- 
lation of infall with time, the source of heavy elements 
in stellar abundances is gradually dominated by those in 
a metal-rich infalling gas. This results in an abundance 
trend directed to the elemental ratios of metal-rich bulge 
stars. Furthermore, an iron contained in a metal-rich 
infalling gas will promote an increase in the [Fe/H] of 
stars, which should result in an extended metal-rich tail 
of ADF. 

4. METAL-RICH OUTFLOW FROM THE GALACTIC 
BULGE 

By incorporation of a switchover of the origin of an in- 
falling gas from the halo to the bulge into the model, we 
describe the evolution of stellar abundances in the Galac- 
tic disk. The basic picture of the model is that the disk 
was formed through an infall of material from outside the 
disk region, and that in the first 9 x lO^yrs, the halo pro- 
vided a very low metallicity gas for the disk; whereas, in 
the subsequent period, a metal-rich gas originating from 
the bulge replaced the former halo gas. The validity of 
such a longtime delay of metal-rich infall should be in- 
vestigated by numerical simulations that will clarify how 
the outflow behaves in the halo under physical processes 
such as a radiative cooling, interactions with halo sub- 
structure, and a possible heating from a prolonged star 
formation in the thick disk (jBensbv et al.l 120071) . For 
simplicity, we propose a formula that is proportional to 
iexp(— t/tin) for the whole duration. In this formula, the 
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Fig. 3. — Features of chemical evolution in the Galactic disk 
predicted by the model in which a metal-rich infall whose elemental 
abundances match those of metal-rich bulge stars are introduced 
in the late stage of disk formation. See the text for details. The 
observational data in the upper panel is same as in Fig. 1. In 
the lower two panel s, th e filled circles and crosse s are taken from 
IBensbv et al.l I I2005I) and lEdvardsson et al.l II1993I '). respectively. 



timecsale and metallicity of an infalling gas is assumed to 
be ti„=5 Gyr and [Fe/H]i„=-1.3 for 0< t(Gyr) <9 and 
ti„=3 Gyr and [Fe/H]i„==+0.2 for 9< t(Gyr) <14. Each 
infall rate is normalized so that an integration over t = 
to oo is an unit. The value of tin for the latter is adjusted 
to give the mass ratio of 4:1 between the former metal- 
poor gas and metal-rich one, which is roughly equivalent 
to the ratio of metal-poor to metal-rich stars in the solar 
neighborhood. We adopt a somewhat smaller [Fe/H]in 
for the bulge origin, taking into account the dilution by 
low-metallicity gas during its passage through the halo. 
Abundances for other heavy elements follow those of ob- 
served halo stars or metal-rich bulge stars for each infall 
(see below). 

Here, we adopt the usual procedure of solving the ba- 
sic equations for chemical evolution, in which the SFR 
and the stellar initial mass function (IMF) are given as 
input parameters. The SFR is assumed to be propor- 



tional to fg with a constant rate coefficient of 0.4 Gyr^ . 
For the IMF, we assume a power-law mass spectrum 
with a Salpeter slope of 1.35, which is combined with 
the nucleosynthesis yie lds of SNe la and II taken from 
iTsuiimoto et al.l (|1995f) . For the hfetime of SN la pro- 
genitors, iia=1.5 Gyr is adopted. 

The results are shown in the three panels of Fig. 3. 
The top panel shows the calculated ADF compared with 
the observations. As shown in the lower two panels, the 
metallicity is enriched up to [Fe/H]'^ +0.3. Despite the 
attainment of super solar metallicity, the predicted ADF 
is not skewed to metal-rich, unlike the models presented 
in Fig. 1. The success in reproducing the right posi- 
tion of the peak, as well as the precise features of both 
sides of its peak, is attributable to the combination of 
a low-metallicty infall with a prolonged timescale at the 
early stage of disk formation, plus the subsequent metal- 
rich infall. The middle panel shows the predicted evo- 
lutionary change in [Mg/Fe] with [Fe/H] for disk stars. 
Here, we assume that the [Mg/Fe] infaii ratio of an in- 
falling gas is -fO.45 for the halo origin and -1-0.3 for 
the bulge origin. First, a delayed supply of iron with 
a negligible amount of Mg from SNe la producs a de- 
creasing trend of [Mg/Fe] for [Fe/H] ^ from a low- 
metallicty infalling gas. Then, it results in [Mg/Fe] ^ at 
[Fe/H]~ for a timescale of about ten billion years. Sub- 
sequently, high abundances of a metal-rich infall, such as 
[Fe/H]infaii=+0.2 and [Mg/H]infaii=+0.5, start to influ- 
ence the stellar abundance patterns so that they continue 
the direction of evolving [Mg/Fe] upward, along with the 
progress of chemical enrichment. 

On the other hand, the lower panel shows that the 
behavior of the [0/Fe] ratio is different. The [0/Fe]infaii 
ratio of metal-rich infall is expected to be nearly the solar 
ratio, as exhibited by metal-rich bulge stars (left middle 
panel of Fig. 2). As a result, such a lower [0/Fe]infaii 
prevents the [0/Fe] trend from turning upward, and in- 
stead it leads to a continuation of the decrease in [O /Fe] 
for [Fe/H] ^ 0, differing from the common upturn fea- 
ture observed with other a-elements. A good agree- 
ment with the observation is achieved by the adoption 
of [0/Fe]infaii=— 0.2, which is somewhat lower than the 
expected ratio. Thus, we predict that larger data sets 
should find the most metal-rich bulge star have lower 
[0/Fe] such as 0.2. 

Finally, we will check whether the supply of metal-rich 
gas satisfies the amount of metal-rick disk stars. Since 
we have neither information on how the metal-rich gas is 
distributed over the disk nor the ADF in the disk other 
than the solar neighborhood, the disk should be con- 
sidered to be one-zone represented by the solar circle. 
Let's take the mass of a bulge as 25% of the mass Mq 
of the Galaxy. The mass fraction n of gas after the star 
formation in the bulge can be deduced from the equa- 
tion Z — phi^~^ (Z: metallicity, p: the effective yield) 
as long as a simple closed-box model is applied. Given 
p deri ved from chemical evo lution in the solar neighbor- 
hood l|Tsuiimoto et al.iri997[ ) , the metallicity of the bulge 
implies /x ~ 0.1. This result combined with the reduc- 
tion in the metallicity of outflow by ~ 0.2 dex through 
mixing with metal-poor halo gas results in ~ 0.05AfG 
for the mass of metal-rich infalling gas. On the other 
hand, the mass 0.75Mg of disk (gas and star) and the 
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fraction 20% of metal-rich ([Fe/H]> 0) stars give the 
mass ^ 0.12Mg of metal-rich stars. Taking into account 
that these stars are made of a metal-rich infall mixed 
with the gas already present in the disk at i = 9 Gyr and 
these contributions are nearly equal in total, the mass of 
metal-rich infalling gas required to make metal-rich disk 
stars is roughly equivalent to the former estimate. 

5. CONCLUSIONS AND DISCUSSION 

We propose that the observational facts of (i) elemen- 
tal abundance features of metal-rich disk stars character- 
ized by an upturn seen for some elements, and (it) the 
extended metal-rich tail up to super solar metallicity in 
the abundance distribution function of disk stars are ev- 
idence that an infall outside the disk region is dominated 
by the gas ejected by the bulge for the recent period of 
five billion years. A grand scheme for the formation of 
bulge and disk in the Galaxy based on this context is de- 
scribed as follows. A prolonged star formation proceeds 
in the bulge, as implied by a clear signature of SNe la 
in the abundances of bulge stars. Then, feedback from 
the supermassive black hole that resides in the Galactic 
center terminates star formation in the bulge and trig- 
gers a highly-enriched outflow. Owing to the relatively 
small mass of the black hole, this outflow is not energetic 
enough to escape from the gravitational potential of the 
Galaxy. Instead, it spends several billion years in the 
halo and then cools down radiatively. Simultaneously, 
the disk is gradually formed through an accretion of ma- 
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angular momentum, timescale, and the mass distribution 
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The discussion of whether it is possible to detect a 
metal-rich infall cannnot be furthered because predict- 
ing its present rate is beyond the capability of our sim- 
plified model. Even so, if we search hydrogen clouds 
already detected in the halo for the potential candidate 
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other hand, high-velocity clouds, such a s Complex C 
types, are likely to be of low metallicity (jWakker et al.l 
iToMlCollins et al.l[2007ll . 
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